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(57) Abstract 

Calculating the mass of air actually flowing into the cylinder with the aid of an inlet pipe filling model which, from the factor angle 
of throttle valve aperture, environmental pressure and parameters representing the valve control, provides a load value on the basis of which 
the injection time is determined. In addition, this load value is used in predictions for estimating the load value at a time at least one 
sampling step later than the actual calculation of the injection time. 



(57) Zusammenfassung 

Berechnung der tatsachlich in den Zylinder einstromenden Luftmasse mit Hilfe eines SaugrohrftUlungsmodells, das aus den 
EingangsgroBen Drosselklappenofrnungswinkel, Umgebungsdruck und Parameter^ die die Ventilsteuerung reprasentieren, cine LastgroBe 
Hefeit, auf dessen Grundlage die Einspritzzeit bestimmt wild. AuSeniem wild diese LastgroBe zur Pradiktion herangezogcn, um die 
LastgroBe zu einem Zeitpunkt abzuschatzen, der mindestens einen Abtastschritt spater liegt als die aktuelle Berechnung der Einspritzzeit. 
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Beschreibung 

Verfahren zum modellgestutzten Bestitnmen der in die Zylinder 
einer Brennkraf tmaschine einstrdmenden Luftmasse 

5 

Die Erfindung betrifft ein Verfahren zum modellgestutzten Be- 
stimmen der in die Zylinder einer Brennkraf tmaschine einstro- 
menden Luftmasse nach dem Oberbegriff des Patentanspruches 1. 

10 Motorsteuerungssysteme fur Brennkraf tmaschinen, die mit 

Kraftstoffeinspritzung arbeiten f benotigen die vom Motor an- 

gesaugte Luftmasse mz y i als ein MaS fur die Motorlast. Diese 
Grofie bildet die Basis zur Realisierung eines geforderten 
Kraftstoff-Luft-Verhaltnisses. Wachsende Anf orderungen an Mo- 

15 torsteuerungssysteme, wie die Verringerung der Schadstoff- 

emission von Kraf tf ahrzeugen, bedingen, dafi die Lastgr&Se fur 
station&re und instationare Vorgange mit geringen zulassigen 
Fehlern bestimmt werden muE. Neben den genannten Betriebsf al- 
ien bietet die genaue Lasterf assung wahrend der Warmlaufphase 

2 0 der Brennkraf tmaschine ein erhebliches Potential zur Schad- 
stof f reduktion . 

Bei luf tmassengefuhrten Motorsteuerungssystemen stellt im In- 
stationarbetrieb das als Lastsignal der Brennkraf tmaschine 

25 dienende Signal des Luf tmassenmessers , der stromauf warts des 
Saugrohrs angeordnet ist, kein MaS fur die tatsachliche Ful- 
lung der Zylinder dar, weil das Volumen des Saugrohrs strom- 
abwarts der Drpsselklappe als Luftspeicher wirkt, der befullt 
und entleert werden muS. Die maSgebende Luftmasse fur die 

30 Einspritzzeitberechnung ist aber diejenige Luftmasse, die aus 
dem Saugrohr heraus und in den jeweiligen Zylinder hinein- 
stromt . 



35 



Bei saugrohrdruckgefuhrten Motorsteuerungssystemen gibt zwar 
das Ausgangssignal des Drucksensors die tatsachlichen Druck- 
verhaltnisse im Saugrohr wieder, die MeSgroSen stehen aber 
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u.a. aufgrund der notwendigen Mittelung der MeSgr&Se erst re- 
lativ spat zur Verfugung. 

Mit der Einfuhrung variabler Ansaugsysteme und variabler Ven- 
txlsteuerungen entstehen fur empirisch gewonnene Modelle zur 
Gewinnung der Lastgrofie aus MeSsignalen eine sehr groEe Viel- 
zahl von EinfluSgrSSen, die die entsprechenden Modellparame- 
ter beeinf lussen. 

Auf physikalischen Ansatzen basierende modellgestutzte Be- 
rechnungsmethoden stellen einen guten Ausgangspunkt zur ge- 
nauen Bestimmung der Luftmasse m Zy i dar. 

Aus der DE 39 19 488 C2 ist eine Vorrichtung zur Regelung und 
zur Vorausbestimmung der Ansaugluf tmenge einer saugrohrdruck- 
gefuhrten Brennkraf tmaschine bekannt, bei der der Drossel- 
klappenof fnungsgrad und die Motordrehzahl als Grundlage zur 
Berechnung des derzeitigen Wertes der in den Brennrautn der 
Maschine eingesaugten Luft verwendet werden. Diese berech- 
nete, gegenwartige Ansaugluf tmenge wird dann als Grundlage 
zur Berechnung des vorausbestimmten Wertes fur die An- 
saugluf tmenge, die in den Brennraum der Maschine zu einer be- 
stimmten Zeit von dem Punkt an, an dem die Berechnung ausge- 
fuhrt wurde, einzusaugen ist, benutzt. Das Drucksignal, das 
stromabwarts der Drosselklappe gemessen wird, wird mitHilfe 
von theoretischen Beziehungen korrigiert, so daE eine Verbes- 
serung der Bestimmung der angesaugten Luftmasse erreicht und 
damit eine genauere Berechnung der Einspritzzeit moglich ist. 

Im instationaren Betrieb der Brennkraf tmaschine ist es aber 
wunschenswert, die Bestimmung der in die Zylinder einstromen- 
den Luftmasse noch genauer durchzufuhren. 

Der Erfindung liegt die Aufgabe zugrunde ein Verfahren anzu- 
geben, mit dem die tatsachlich in den Zylinder der Brenn- 
kraf tmaschine einstrdmende Luftmasse mit hoher Genauigkeit 
bestimmt werden kann. AulSerdem sollen systembedingte Totzei- 
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ten, die aufgrund der Kraf tstof fvorlagerung und der Rechen- 
zeit bei der Berechnung der Einspritzzeit auftreten konnen, 
kompensiert werden, 

5 Diese Aufgabe wird gemaS den Merkmalen des Patentanspruches 1 
gelost . 

Vorteilhafte Weiterbildungen finden sich in den Unteranspru- 
chen . 

10 

Ausgehend von einem bekannten Ansatz ergibt sich eine Modell- 
beschreibung, die auf einer nichtlinearen Differential - 
gleichung basiert. Im folgenden wird eine Approximation die- 
ser nichtlinearen Gleichung vorgestellt. Ira Ergebnis dieser 

15 Approximation lafit sich das Systemverhalten mittels einer bi- 
linearen Gleichung beschreiben, die die schnelle Losung der 
Beziehung im MotorsteuergerSt des Kraf tf ahrzeugs unter Echt- 
zeitbedingungen gestattet. Der gewahlte Modellansatz beinhal- 
tet dabei die Modellierung von variablen Saugsystemen und Sy- 

20 stemen mit variablen Ventilsteuerungen. Die durch diese An- 

ordnung und durch dynamische Nachladung, d.h. durch Reflexio- 
nen von Druckwellen im Saugrohr hervorgeruf enen Effekte, kon- 
nen ausschliefilich durch die Wahl stationar bestimmbarer Pa- 
rameter des Modelles sehr gut berucksichtigt werden. Alle Mo- 

25 dellparameter sind einerseits physikalisch interpretierbar 

und andererseits ausschliefilich aus station^ren Messungen zu 
gewinnen , 

Die meisten Algorithmen zur zeitdiskreten Losung der Diffe- 
30 rentialgleichung, die das Verhalten des hier genutzten Model- 
les beschreibt, erfordern vor allem bei geringem Druckabfall 
uber der Drosselklappe, d.h. bei Vollast eine sehr kleine Re- 
chenschrittweite, um numerisch stabil zu arbeiten. Die Folge 
ware ein unvertretbarer Rechenauf wand bei der Bestimmung der 
35 LastgroSe. Da Lasterf assungssysteme meist segment synchron ar- 
beiten, d.h. fur 4-Zylindermotoren wird alle 180° KW ein MeS- 
wert abgetastet, muS die Modellgleichung ebenfalls segment- 
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synchron gelost werden. Im nachfolgenden wird ein absolut 
stabiles Differenzenschema zur Losung von Dif ferentialglei- 
chungen eingesetzt, das numerische Stability bei beliebiger 
Schrittweite garantiert. 



Das erfindungsgemafie modellgestiitzte Berechnungsverf ahren 
bietet zudem die Moglichkeit einer Pradiktion des Lastsigna- 
les urn eine wahlbare Anzahl von Abtastschritten, d.h. eine 
Vorhersage des Lastsignales mit variablem Pradiktionshori- 
zont. Wird die dem Pradiktionshorizont bei konstanter Dreh- 
zahl proportionale Pradiktionszeit nicht zu grofi, so erhalt 
man ein pradiziertes Lastsignal hoher Genauigkeit. 

Eine solche Vorhersage ist notwendig, da zwischen der Erfas- 
sung relevanter Mefiwerte und der Berechnung der LastgrofSe ei- 
ne Totzeit entsteht. Desweiteren mufi aus Grunden der Ge- 
mischaufbereitung vor dem eigentlichen Beginn der Ansaugphase 
des jeweiligen Zylinders moglichst genau die Kraf tstof fmasse 
iiber die Einspritzventile zugemessen werden, die im Verlauf 
der kommenden Ansaugphase im gewunschten Verhaltnis zur Luft- 
masse m Zyl steht. Ein variabler Pradiktionshorizont verbessert 
die Gute der Kraf tstof fzumessung im instationaren Motorbe- 
trieb. Da bei steigender Drehzahl die Segmentzeit abnimmt, 
muf§ der Einspritzvorgang eine groSere Anzahl von Segmented 
eher beginnen, als dies bei einer niedrigeren Drehzahl der 
Fall ist. Urn die zu dosierende Kraf tstof fmasse moglichst ex- 
akt bestimmen zu kdnnen, ist die Pradiktion der Lastgrofie urn 
die Anzahl von Segmenten, urn die die Kraf tstof fvorlagerung 
vorgenommen wird, notwendig, urn ein gefordertes Kraftstoff- 
Luft-Verhaltnis auch in diesem Fall einzuhalten. Die Pradik- 
tion der Lastgrofie tragt somit aus einer wesentlichen Verbes- 
serung der Einhaltung des gefordertene Kraf tstof f-Luft- Ver- 
haltnisses im instationaren Motorbetrieb bei. Dieses System 
zur modellgestutzten Lasterf assung ist in den bekannten Mo- 
torsteuerungssystemen, d.h. bei luf tmassengefiihrte bzw. saug- 
rohrdruckgefuhrte Motorsteuerungssysteme wird im folgenden 
ein Korrekturalgorithmus in Form eines Modellregelkreises 
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formuliert, der bei auf tretenden Ungenauigkeiten von Modell- 
parame tern eine permanente Genauigkeitsverbesserung, d.h. ei- 
nen Modellabgleich im stationSren und instationaren Betrieb 
gestattet . 

Ein Ausfiihrungsbeispiel des erf indungsgemaEen Verfahrens wird 
anhand der nachf olgenden schematischen Zeichnungen beschrie- 
ben. Dabei zeigen: 



Figur 1 eine Prinzipskizze zum Saugsystem einer Otto-Brenn- 
kraf tmaschine einschlieSlich der entsprechenden Mo- 
dell- und MefigroSen, 

Figur 2 die Durchf lufif unktion und die dazugehorige Polygon- 
zugapproximation, 

Figur 3 eine Prinzipdarstellung zum Modellregelkreis fur 
luf tmassengefuhrte Motors teuerungssyst erne und 

Figur 4 eine Prinzipdarstellung zum Modellregelkreis fur 
saugrohrdruckgef uhrte Motorsteuerungssysteme . 

Bei der modellgestiitzten Berechnung der Lastgrofie mz y i wird 
von der in Figur 1 dargestellten prinzipiellen Anordnung aus- 
gegangen. Aus Griinden der Ubersichtlichkeit wird dabei nur 
ein Zylinder der Brennkraf tmaschine dargestellt. Mit dem Be- 
zugzeichen 10 ist dabei ein Saugrohr einer Brennkraf tmaschine 
bezeichnet, in dem eine Drosselklappe 11 angeordnet ist. Die 
Drosselklappe 11 ist mit einem, den Offnungsgrad der Drossel- 
klappe ermittelnden Drosselklappenstellungsf uhler 14 verbun- 
den. Stromauf warts der Drosselklappe 11 ist bei einem luft- 
massengefuhrten Motorsteuerungssystem ein Luf tmassenmesser 12 
angeordnet, wahrend bei einem saugrohrdruckgef uhrten Motor- 
steuerungssystem ein Saugrohrdruckfuhler 13 im Saugrohr ange- 
ordnet ist. Je nach Art der Lasterf assung ist somit nur eine 
der beiden Komponenten 12, 13 vorhanden. Die Ausgange des 
Luf tmassenmessers 12, des Drosselklappenstellungsgebers 14 
und des zum Luf tmassenmesser 12 alternativ vorhandenen Saug- 
rohrdrucksensors 13 sind mit Eing&ngen einer nicht darge- 
stellten, an sich bekannten elektronischen Steuerungseinrich- 
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tung der Brennkraftmaschine verbunden. AuSerdem sind in Figur 
1 noch schematisch ein EinlaSventil 15, ein AuslaSventil 16, 
sowie ein in einem Zylinder 17 beweglichen Kolben 18 darqe- ' 
stellt. 

AuSerdem sind in Figur 1 ausgewahlte GroSen bzw. Parameter 
des Saugsystems eingezeichnet . Dabei bedeutet das Dachsymbol 
"a» tiber einer GroSe, daS es sich urn eine ModellgroSe han- 
delt, wahrend GrSSen ohne Dachsymbol "a" MeSgroSen reprasen- 
tieren. Im einzelnen bedeuten: 

PU Umgebungsdruck, P s Saugrohrdruck , T s Temperatur der Luft 
im Saugrohr, V s das Volumen des Saugrohrs. 

GrdSen mit einem Punktsymbol kennzeichnen die erste zeitliche 
Ableitung der entsprechenden Gr&Sen. m DK ist somit der Luft- 
massenstrom an der Drosselklappe und m Zy , ist der Luftmassen- 
strom der tatsachlich in den Zylinder der Brennkraftmaschine 
einstromt. 

Die grundlegende Aufgabe bei der mcdellgestutzten Berechnung 
des Motorlastzustandes besteht nun in der Losung der Diffe- 
rent ialgleichung fur den Saugrohrdruck 



D T ( A A \ 

m DK~ m Zyl 



V 5 



(2.1) 



die sich unter der Voraussetzung konstanter Temperatur der 
Luft im Saugrohr T s aus der Zustandsgleichung idealer Case 
herleiten laSt . 

Mit R L ist dabei die allgemeine Gaskonstante bezeichnet. 

A 

Die Lastgrofie m Zy l wird durch Integration aus dem Zylinder- 

A 

massenstrom m Zy i bestimmt. Die durch (2.1) beschriebenen 
Verhaltnisse sind auf Mehrzylinder-Brennkraf tmaschinen mit 
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Schwingrohr- (Schaltsaugrohr- ) und/oder Resonanzsaugsysteme 
ohne strukturelle Anderungen anwendbar. 

Fur Systeme mit Multi-Point-Einspritzungen, bei denen die 
Kraf tstof f zumessung durch mehrere Einspritzventile erfolgt, 
gibt die Gleichung (2.1) die Verhaltnisse genauer wieder als 
dies bei Single-Point -Einspritzungen, d.h. bei Einspritzun- 
gen, bei denen der Kraftstoff mittels eines einzigen Kraft- 
stof feinspritzventiles zugemessen wird, der Fall ist. Bei 
erstgenannter Art der Kraf tstof f zumessung ist nahezu das ge- 
samte Ansaugsystem mit Luft gefullt. Lediglich in einem klei- 
nen Bereich vor den EinlaSventilen befindet sich ein Kraft- 
stoff -Luftgemisch . Im Gegensatz dazu ist bei Single-Point - 
Einspritzsystemen das gesamte Saugrohr von der Drosselklappe 
bis zum EinlaSventil mit Kraf tstof f -Luf t-Gemisch gefullt, da 
das Einspritzventil vor der Drosselklappe angeordnet ist. In 
diesem Fall stellt die Annahme eines idealen Gases eine star- 
kere Naherung dar, als dies bei der Multi-Point-Einspritzung 
der Fall ist. Bei Single-Point -Einspritzung erfolgt die 
Kraf tstof f zumessung entsprechend 
? a 
moK , bei Multi-Point-Einspritzung entsprechend mzyt . 

A 

Im folgenden wird die Berechnung der Massenstrome mDK und 

A 

mzyi naher beschrieben. 

Die ModellgrolSe des Luf tmassenstromes an der Drosselklappe 

A 

moK wird durch die Durchf lufigleichung idealer Gase durch 
Drosselstellen beschrieben. An der Drosselstelle auftretende 
Stromungsverluste werden durch den reduzierten Stromungsquer- 

A A 

schnitt Ared beriicksichtigt . Der Luf tmassenstrom moK wird 
demnach durch die Beziehung 
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mit 



2 



\ p v ) 



fur uberkritische Druckverhaltnisse 



bzw. 



V * const, fur kritische Druckverhaltnisse 



bestimmt 



(2.2) 



A 

A RED 
K : 

R L : 

A 

Pu: 

A 

V : 



ModellgroSe des Luf tmassenstromes an der Drossel- 
klappe 

reduzierter Stromungsquerschnitt 
Adiabatenexponent 
allgemeine Gaskonstante 
Temperatur der Luft im Saugrohr 

ModellgroSe des Umgebungsdruckes 

Modellgrofie des Saugrohrdruckes 
Durchf lufif unkt ion . 



An der Drosselstelle, d.h. an der Drosselklappe auftretende 
Stromungsverluste werden uber die geeignete Wahl von Ared 
berucksichtigt, Aus stationaren Messungen kann bei bekannten 
Drucken vor und hinter der Drosselstelle und bekanntem Mas- 
senstrom durch die Drosselstelle eine Zuordnung zwischen dem 
vom Drosselklappenstellungsfuhler 14 ermittelten Drosselklap- 
penwinkel und dem entsprechendem reduzierten Querschnitt 

A 

Ared angegeben werden. 



A 

Wird der Luf tmassenstrom ntDK an der Drosselklappe durch die 
Beziehung (2.2) beschrieben, so entsteht ein kotnplizierter 
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Algorithmus zur numerisch richtigen Losung der Dif f erential- 
gleichung (2.1). Zur Reduktion des Rechenaufwandes wird die 
Durchf lufcfunktion \\t durch einen Polygonzug approximiert . 

5 Figur 2 zeigt den Verlauf der Durchf luSfunktion \y und das 

darauf angewandte Approximationsprinzip. Innerhalb eines Ab- 
schnittes i (i = l...k) wird die Durchf luSfunt ion y durch 
eine Gerade dargestellt. Mit einer vertretbaren Anzahl von 
Geradenabschnitten kann damit eine gute Approximation er- 
10 reicht werden. Durch einen solchen Ansatz kann die Gleichung 
(2.2) zur Berechnung des Massenstromes an der Drosselklappe 

A 

mDK durch die Beziehung 



15 



/ 2k 1 



moK^APPRox = Ared-J-t- Pu 



Ps 



(2.3) 



fur i = (l. . .k) 
approximiert werden. 

In dieser Form beschreibt mi die Steigung und ni das Absolut- 
20 glied des jeweiligen Geradenabschnittes . Die Werte fur die 
Steigung und fur das Absolutglied werden in Tabellen als 
Funktion des Verhaltnisses Saugrohrdruck zu Umgebungsdruck 

A 

abgelegt. 

25 Auf der Abszisse von Figur 2 ist dabei das Druckverhaltnis 

A 

•— und auf der Ordinate der Funktionswert {0 - 0.3) der 
Pu 

Durchf luSfunkt ion y aufgetragen. 

A K 

Fur Druckverhaltnisse <[ -]* ' ist y = konstant, d.h. 

u 

30 daS der Durchf luS an der Drosselstelle nur noch vom Quer- 

schnitt abhangig ist und nicht mehr von den Druckverhaltnis- 
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sen Die in die jeweiligen Zylinder der Brennkraf tmaschine 
einstromende Luftmasse lafit sich analytisch nur schwer be- 

de'rTunL 8 " " ^s^ 1 ^gt . Die Fullung 

der Zylinder W1 rd wextestgehend durch den Saugrohrdruck di/ 
Drehzahl und durch die Ventilsteuerzeiten bestimmt 

Zur moglichst genauen Berechnung des Massenstroms in den je- 
weiligen Zylinder l Zyl ist deshalb einerseits die Beschrei- 
bung der VerMltnisse i m Ansaugtrakt der Brennkraf tmaschine 
1 B P h 3rtleller ^"erentialgleichungen und andererseits 
die Berechnung des Massenstromes am EinlaSventil nach der 
DurchfluSgleichung als erf orderliche Randh*^ 

Er8 . . . ,7 . UCfie Randbedingung notwendig. 

Erst dieser komplizierte Ansatz gestattet die Berucksichti- 
gung dynamischer Nachladeef fekte, die von der Drehzahl der 
Saugrohrgeometrie, der Zylinderzahl sowie den Ventilsteuer- 
zeiten maSgeblich beeinfluSt werden. 

Da eine Berechnung nach oben genanntem Ansatz in der elektro 
nischen Steuerungseinrichtung der Brennkraf tmaschine nicht 
realisierbar ist, geht eine mogliche Naherung von einem ein- 
fachen Zusammenhang zwischen Saugrohrdruck P s und Zyindermas- 
senstrom m Zyl aus . Fur einen weiten Bereich der sinnvollen 
Ventilsteuerzeiten kann dafur in guter Naherung von einem li- 
nearen Ansatz der Form 

A 

A 

mzyi_ A ppRox =y r P s + y Q ^ 
ausgegangen werden. 

Die Steigung y, und das Absolutglied y 0 der Beziehung (2 4) 
Bind dabei, unter Berucksichtigung aller wesentlichen Ein- 
fluSfaktoren Funktionen der Drehzahl, der Saugrohrgeometrie 
der Zylinderzahl, der Ventilsteuerzeiten sowie der Temperatur 
der Luft im Saugrohr T s . Die Abhangigkeit der Werte von y, 
und y 0 von den Einf luSgrSSen Drehzahl, Saugrohrgeometrie, Zy- 
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linderzahl und den Ventilsteuerzeiten und Ventilerhebungskur- 
ven kann dabei uber stationare Messungen ermittelt werden. 
Uber diese Wertebestimmung wird ebenfalls der EinluS von 
Schwingrohr- und/oder Resonanzsaugsystemen auf die von der 
5 Brennkraf tmaschine angesaugte Luftmasse gut wiedergegeben. 
Die Werte von y, und y 0 sind in Kennfeldern der elektroni- 
schen Motorsteuerungseinrichtung abgelegt. 

Als bestimmende GroSe zur Ermittlung der Motorlast wird der 
10 Saugrohrdruck P s ausgewahlt. Mit Hilfe der Modell-Dif f eren- 
tialgleichung soil diese Grdfie moglichst exakt und schnell 

A 

geschatzt werden. Die Schatzung von Ps erfordert die Losung 
der Gleichung (2.1). 



15 



Mit den anhand der Formeln (2.2) und (2.3) eingefuhrten Ver- 
einfachungen kann (2.1) durch die Beziehung 



20 



25 



A 

A RED 



I 2k 1 



A 

■Pu- 



A 

m i • +w / 

Pu 



-jri Ps +y 0 J 



fur i = (1. . .k) 



(2.5) 



approximiert werden. Betrachtet man, entsprechend den Voraus- 
setzungen zur Herleitung von Gleichung (2.1), die Temperatur 
der Luft im Saugrohr T s als eine langsam veranderliche Mefi- 

A 

groSe sowie Ared als EingangsgrdSe, so lafit sich die nicht- 
lineare Form der Dif f erentialgleichung (2.1) durch die bili- 
neare Gleichung (2.5) approximieren . 



30 



Zur Losung der Gleichung (2.5) wird diese Beziehung in eine 
geeignete Dif f erenzengleichung ubergefuhrt. 



35 



Als Kriterium zur Auswahl des geeigneten Dif f erenzenschemas 
konnen die folgenden prinzipiellen Anf orderungen an die L6- 
sungseigenschaf ten der zur bildenden Dif f erenzengleichung 
formuliert werden: 



WO 96/32579 



PCT/DE96/00615 



12 



1. Das Differenzenschema muE auch unter extremen dynamischen 
Anforderungen konservativ sein, d.h. die Losung der Diffe- 
renzengleichnung mu£ der Losung der Dif f erentialgleichung 
entsprechen, 

2. die numerische Stabilitat muS zu Abtastzeiten, die den ma- 
ximal moglichen Segmentzeiten entsprechen, im gesamten Ar- 
beit sbereich das Saugrohrdruckes garantiert sein. 

Forderung 1 ist durch einen impliziten Rechenalgorithmus er- 
fiillbar. Aufgrund der Approximation der nichtlinearen Diffe- 
rentialgleichung (2.1) durch eine bilineare Gleichung ist das 
entstehende implizite Ldsungsschema ohne Einsatz iterativer 
Verfahren losbar, da die Dif f erenzengleichung in eine expli- 
zite Form dberfiihrt werden kann. 

Die zweite Forderung ist aufgrund der Konditionierung der 
Dif f erentialgleichung (2.1) und deren Approximation (2.5) nur 
durch eine Rechenvorschrif t zur Bildung der Differenzen- 
gleichung erfiillbar, die absolut stabil arbeitet. Diese Ver- 
fahren werden auch als A-stabile Verfahren bezeichnet. Kenn- 
zeichnend fur diese A- Stabilitat ist die Eigenschaft des Al- 
gorithms, bei einem stabilen Ausgangsproblem fur beliebige 
Werte der Abtastzeit, d.h. Segmentzeit T A numerisch stabil zu 
sein. Eine mogliche Rechenvorschrif t zur numerischen Losung 
von Differentialgleichungen, die beiden Forderungen gerecht 
wird, ist die Trapezregel . 

Die durch Anwendung der Trapezregel entstehende Diffe- 
renzengleichung lautet im vorliegenden Fall 

P s [N) = P S [N-\] + Z~U [JV-lJ + P* [N] 

V 

fur N = il..cP) (2.6) 
def iniert . 
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Wird diese Vorschrift auf (2.5) angewandt, so ergibt sich die 
Beziehung 



Ps [»] = 



Ps[*-lh%~Ps [N-\] 



\-I±.ELlIk. 

2 V, 



| 2K 1 



T s 



T A R L T S 



*RED 



2k 1 



k-1 R L -T S 



Pu"i -Yo 



2 f 5 ^ ^ 



10 fiir N = und i = (l...k) (2.7) 

zur Berechnung des Saugrohrdruckes P s [N] als MaS fur die Mo- 
torlast . 



[N] bedeutet dabei das aktuelle Segment bzw. der aktuelle Re- 
chenschritt, + das nachstf olgende Segment bzw. der 
nachstf olgende Rechenschritt . 



Im folgenden wird die Berechnung des aktuellen und pradizier- 
ten Lastsignales beschrieben. 

20 

A 

Aus dem berechneten Saugrohrdruck Ps kann der Luftmassen- 

A 

Strom mzyi der in die Zylinder einstromt, durch die Beziehung 
(2.4) ermittelt werden. Wendet man einen einfachen Integra - 
tionsalgorithmus an, so erhalt man fur die wahrend eines An- 
25 saugtaktes von der Brennkraf tmaschine angesaugte Luftmasse 
die Beziehung 

a 7* ( A A ^ 

m Zy i[N) = -±-\ mzyi [N-\] + m Zy i [N] 



30 



fur N = (l. . Xfi) 



(2.8) 
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Dabei wird davon ausgegangen, daS der Anfangswert der Last- 
groSe null ist . Fur die segment synchrone Lasterf assung sinkt 
mit steigender Drehzahl die Segment zeit , wahrend die Segment - 
anzahl, urn die eine Kraf tstof fvorlagerung vorgenommen wird, 
5 steigen muS. Aus diesem Grund ist es erf orderlich, die Pra- 

diktion des Lastsignals fur einen verSnderlichen Pradiktions- 
horizont H, d.h. fur eine bestimmte, in erster Linie dreh- 
zahlabhangige Anzahl H von Segmenten, auszulegen. Berucksich- 
tigt man diesen verander lichen Pradiktionshorizont H, so kann 
10 Gleichung (2.8) in der Form 

a T ( * A ^ 

m Z yi[N + H] = -±- mzyi [AT + //-1]-h m Zy t [N + H] 

2 V ) 

fur N = (1. . <fi) (2.9) 

15 geschrieben werden. 

Fur die weiteren Uberlegungen wird davon ausgegangen, daE 
sich die Segmentzeit T A und die Parameter y, und y 0 der Be- 

A 

ziehung (2.4), die zur Bestimmung des Massenstromes mz y i aus 

A 

20 dem Saugrohrdruck P s erforderlich sind, liber die Pradiktions- 
zeit nicht andern. 

Unter dieser Voraussetzung wird die Pradiktion eines Wertes 

A 

fur mz y i[N + H] durch die Pradiktion des entsprechenden Druck- 

A 

wertes Ps [N + H] erreicht. Dadurch nimmt die Gleichung (2.9) 
25 die Form 



mzyi[ 

fur N = (1. . .<X>) (2.10) 
3 0 an. 

Da bei dem beschriebenen Verfahren die zeitliche Anderung des 

A 

Saugrohrdruckes Ps in analytischer Form vorliegt, wird im 

A 

folgenden die Pradiktion des Druckwertes Ps [N + H] durch H- 
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fache Anwendung der Trapezregel erreicht. In diesem Fall er- 
halt man die Beziehung 



Ps [N + H]= P S [N] + ^-H 



Ps[*-1]+P,[ 



fttr N = (i. .00 ) 



(2.11) 



10 



A 

Bestimmt man den Druck Ps [AT + ff-l] in analoger Weise, so kann 
fiir das pradizierte Lastsignal die Gleichung 



fur N « <1. .GO) 



(2.12) 



15 



angegeben werden. 

Wahlt man fiir den Pradiktionshorizont H Werte in der Grofien- 
ordnung von 1...3 Segmenten, so kann .mit der Formel (2.12) 
ein gut pradiziertes Lastsignal erhalten werden. 



20 Im folgenden wird das Prinzip des Modellabgleichs fiir luft- 
massen- und saugrohrdruckgef uhrte Motorsteuerungssysteme er- 
klart . 



25 



30 



Bedingt durch den Einsatz von Motoren mit variabler Ventil- 
steuerung und/oder veranderlicher Saugrohrgeometrie, durch 
Fertigungstoleranzen und Alterungserscheinungen, sowie durch 
Temperatureinf lusse sind die Werte von y } und y 0 mit einer 
gewissen Unsicherheit behaftet. Die Parameter der Gleichung 
zur Bestimmung des Massenstromes in den Zylindern sind, wie 
oben beschrieben, Funktionen vielfal tiger Einf luSgroEen, von 
denen nur die wichtigsten erfaSt werden konnen. 



35 



Bei der Berechnung des Massenstromes an der Drosselklappe 
wirken sich MeSfehler bei der Erfassung des Drosselklappen- 
winkels und Approximationsf ehler bei der Polygonzugapproxima- 
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tion der Durchf luSf unktion y auf die Modellgrofien aus. Be- 
sonders bei kleinen Drosselklappenwinkeln ist die System- 
etnpf indlichkeit gegenuber erstgenannten Fehlern besonders 
hoch. Daraus ergibt sich, daS kleine Anderungen der Drossel- 
5 klappenstellung einen gravierenden EinluS auf Massenstrom 

bzw. Saugrohrdruck haben. Urn die Wirkung dieser Einflusse zu 
reduzieren, wird im folgenden ein Verfahren vorgeschlagen, 
das es gestattet, bestimmte Grofien, die EinfluS auf die Mo- 
dellrechnung haben, so zu korrigieren, daS eine genauigkeits 
10 verbessernde Modellanpassung fur stationaren und instationa- 
ren Motorbetrieb durchgefuhrt werden kann. 

Die Anpassung wesentlicher Parameter des Modells zur Bestim- 
mung der LastgroSe der Brennkraf tmaschine erfolgt durch die 
15 Korrektur des aus dem gemessenen Drosselklappenwinkel be- 

A 

stimmten reduzierten Querschnitts Ared durch die Korrektur- 

A 

groSe A A red . 

Die Eingangsgrofie zur korrigierten Saugrohrdruckberechnung 

A 

20 Ared wird damit durch die Beziehung 



A A A 

A redkorr = Ared + A Ared (3.11) 
25 beschrieben. 

In der Gleichung (2.2) und nachf olgenden Formeln wird dann 

A A 

Ared durch Aredkorr ersetzt. Zur Verbesserung des Folgever- 
haltens des Regelkreises wird der aus dem MeSwert des Dros- 
30 selklappenwinkels abgeleitete reduzierte Drosselklappenquer- 

A 

schnitt Ared in die Model 1 re chnung einbezogen. Die Korrek- 

A 

turgroSe AAred wird durch Realisierung eines Modellregel- 
kreises gebildet. 
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Fur luftmassengefuhrte Motorsteuerungssysteme ist der mittels 
des Luftmassenmessers an der Drosselklappe gemessene Luft- 
massenstrom hidk_lmm die FuhrungsgroSe dieses Regelkreises , 
w&hrend fdr saugrohrdruckgefiihrte Systeme der gemessene Saug- 
rohrdruck P s als FuhrungsgroSe genutzt wird. Uber eine Fol- 

A 

geregelung wird der Wert von AAred so bestimmt, daS die Re- 
gelabweichung zwischen FuhrungsgroSe und der ensprechenden 
Regelgr6Se minimiert wird. 

Urn auch im dynamischen Betrieb Genauigkeitsverbesserungen mit 
der genannten Methode zu erreichen, muS die MeSwerterf assung 
der FuhrungsgroSe moglichst exakt nachgebildet werden. In den 
meisten Fallen sind dabei das dynamische Verhalten des Sen- 
sors, d.h. entweder des Luftmassenmessers oder des Saugrohr- 
drucksensors und eine nachfolgend durchgef uhrte Mittelwert- 
bildung zu berucksichtigen. 



Das dynamische Verhalten des jeweiligen Sensors kann in er- 
ster Naherung als ein System erster Ordnung mit eventuell ar- 
20 beitspunktabhangigen Verzogerungszeiten T± modelliert werden. 
Im Falle eines luf tmassengefuhrten Systems lautet eine mdg- 
liche Gleichung zur Beschreibung des Sensorverhaltens 



25 



m D K_LMM [N] = e r ' • rriDK + 



-id 



m D K_iMM[N-\] (3.12) 



Eine GroSe, die beim gewahlten Ansatz einen wesentlichen Ein- 

A 

fluS auf den maximal moglichen Massenstrom mz y \ besitzt, ist 

A 

der Umgebungsdruck Pv . Aus diesem Grund kann nicht von einem 
konstanten Wert dieser GroSe ausgegangen werden, sondern es 
3 0 erfolgt eine Anpassung in der nachfolgend beschriebenen Art 
und Weise. 

A 

Der Wert des Umgebungsdruckes Pv wird verSndert, wenn der 
Betrag der KorrekturgroSe AAred eine bestimmte Schwelle 
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uberschreitet oder wenn das Druckverhaltnis groEer als 

i 

eine wahlbare Konstante ist. Damit wird gewahrleistet , daS 
sowohl im Teil- als auch im Vollastbereich eine Umgebungs- 
druckanpassung erfolgen kann. 

5 

Im folgenden wird ein Modellabgleich fur luf tmassengefuhrte 
Motorsteuerungssysteme erklart. Fur dieses System kann die in 
Figur 3 dargestellte Modellstruktur angegeben werden. 

10 Der Drosselklappenstellungsfuhler 14 (Figur 1) liefert ein 

dem Offnungsgrad der Drosselklappe 11 entsprechendes Signal, 
z.B. einen Drosselklappenof f nungswinkel . In einem Kennfeld 
der elektronischen Motorsteuerungseinrichtung sind zu ver- 
schiedenen Werten dieses Drosselklappenof fnungswinkels zuge- 

15 horige Werte fur den reduzierten Querschnitt der Drossel- 

A 

klappe Ared abgespeichert . Diese Zuordnung wird durch den 
Block "statisches Modell" in Figur 3 und in Figur 4 reprasen- 
tiert. Das Teilsystem "Saugrohrmodell 11 in den Figuren 3 und 4 
reprasentiert das durch (2.7) beschriebene Verhalten. Fuh- 

20 rungsgrofie dieses Modellregelkreises ist der MeSwert des iiber 
ein Segment gemittelten Luf tmassenstromes an der Drossel- 
klappe moK^LMM ■ Wird als Regler in diesem Model lregelkreis 
ein PI -Regler eingesetzt, so ist die bleibende Regelab- 
weichung null, d.h. ModellgroSe und MeSgroSe des Luftmassen- 

25 stromes an der Drosselklappe sind identisch. 

Die Pulsationserscheinungen des Luf tmassenstromes an der 
Drosselklappe, die vor allem bei 4 -Zylindermotoren zu be- 
obachten sind, fuhren bei betragsbildenden Luf tmassenmessern 
zu erheblichen positiven MeSfehlern und somit zu einer stark 

30 f ehlerbehaf teten FuhrungsgroSe . Durch eine Abschaltung des 
Reglers, d.h. einer Verkleinerung der Reglerparameter kann 
zum gesteuerten modellgestutzten Betrieb ubergegangen werden. 
Bereiche, in denen die genannten Pulsationen auftreten, kon- 
nen somit mit dem selben Verfahren unter Berucksichtigung dy- 

35 namischer Zusammenhange behandelt werden, wie diejenigen Be- 
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reiche, in denen eine nahezu ungestorte Fuhrungsgr6£e vor- 
liegt. Im Gegensatz zu Verfahren, die relevante MeSwerte nur 
in stationaren Betriebspunkten berucksichtigen, bleibt das 
beschriebene System nahezu uneingeschrankt arbeitsf ahig. Bei 
Ausfall des Luf tmassensignals oder des Signals des Drossel- 
klappenstellungsfuhlers ist das vorgestellte System in der 
Lage, ein entsprechendes Ersatzsignal zu bilden. Bei Ausfall 
der FuhrungsgroSe muS der gesteuerte Betrieb realisiert wer- 
den, wahrend im anderen Fall der geregelte Betrieb die kaum 
beeintrachtigte Funktionsf ahigkeit des Systems garantiert. 

Der Block "Saugrohrmodell" reprasentiert die Verhaltnisse wie 
sie anhand der Gleichung (2.7) beschrieben sind und hat dem- 

zufolge als AusgangsgroGe die Modellgrofie P s sowie die zeit- 

A A 

liche Ableitung P s und die Grofie ttidk . Nach der Modellierung 
des Sensoriibertragungsverhaltens d.h. des Ubertragungsverhal- 
tens des Luf tmassenmessers und der Abtastung wird die Modell- 

groSe ntDK^LMM einer Mittelung unterzogen, so daS die gemit- 

A 

telte GroEe ttidk^imm und der vom Luf tmassenmesser gemessene 

durchschnittliche Luf tmassenstrom moK^LMM einem Vergleicher 
zugefuhrt werden konnen. Die Differenz beider Signale bewirkt 

A 

eine Anderung AAred des reduzierten Strdmungsquer-schnittes 

A 

Ared, so daS stationar und instationar ein Mo-dellabgleich 
erfolgen kann. 

Fur saugrohrdruckgefuhrte Motorsteuerungssysteme wird die in 
Figur 4 dargestellte Modellstruktur angegeben, wobei gleiche 
Blocke wie in Figur 3 gleiche Bezeichnungen tragen. Ebenso 
wie bei dem luf tmassengef uhrten Motorsteuerungssystem repra- 
sentiert das Teilsystem "Saugrohrmodell", das durch die Dif- 
ferenzengleichung {2.1) beschriebene Verhalten. FuhrungsgroSe 
dieses Modellregelkreises ist der MeSwert des uber ein Seg- 
ment gemittelten Saugrohrdruckes Ps_s . Wird ebenfalls wie in 
Figur 3 ein PI-Regler eingesetzt, so ist im stationaren Fall 
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der MeSwert des Druckes im Saugrohr Ps_s mit der ModellgroSe 

A 

Ps_s identisch. Wie oben beschrieben, bleibt auch das vorlie- 
gende System nahezu uneingeschrankt arbeitsf ahig, da bei Aus- 
fall des Saugrohrdrucksignales oder des MeSwertes fur den 
Drosselklappenwinkel ein entsprechendes Ersatzsignal gebildet 
werden kann. 



Die durch das Saugrohrmodell erhaltenen ModellgroSen P s , P s 
werden einem Block "Pradiktion" zugefiihrt. Da mit den Model- 
10 len auch die Druckanderungen im Saugrohr berechnet werden, 
konnen diese Druckanderungen dazu verwendet werden, den zu- 
kunftigen Druckverlauf im Saugrohr und damit die Zylinder- 
luftmasse fur das nachste + oder fur die nachsten Seg- 
ment e [N + H] zu schatzen. Die Grofie mz y i bzw. die GroSe 

15 mzyi[N + \] dienen dann zur exakten Berechnung der Einspritz- 
zeit, wahrend derer Kraftstoff eingespritzt wird. 
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Patent anspruche 

1. Verfahren zum Bestimraen der einstromenden Luftmasse in den 
bzw. die Zylinder einer Brennkraf tmaschine mit 
5 - einem Ansaugsystem, das ein Saugrohr (10) und eine darin 
angeordnete Drosselklappe (11) , sowie einen den Offnungs- 
grad der Drosselklappe (19) erfassenden Drosselklappenstel- 
lungsfuhler (14) aufweist, 

- einem, ein Lastsignal (m D K_iMMlPs_s) der Brennkraf tmaschine 
10 erzeugenden Sensor (12; 13) 

- einer elektrischen Steuerungseinrichtung, die auf der 

Grundlage des gemessenen Lastsignals (mDK_LMM\Ps_s) und der 
Drehzahl der Brennkraf tmaschine eine Grundeinspritzzeit be- 
rechnet , 

15 dadurch gekennzeichnet, daS 

- die Verhaltnisse im Ansaugsystem mittels eines Saugrohrful- 
lungsmodells nachgebildet werden, wobei als EingangsgroSen 
des Modells der Offnungsgrad der Drosselklappe (11) , der 
Umgebungsdruck (P v ) und die Ventilstellung reprasentieren- 

20 de Parameter herangezogen werden, 

A 

- eine ModellgroSe fur den Luf tmassenstrom (^)an der Dros- 
selklappe (11) mit Hilfe der Durchf luSgleichung idealer Ga- 
se durch Drosselstellen beschrieben wird (Gleichung 2.2) , 

A 

- eine ModellgroSe fur den Luf tmassenstrom (mz y i) in den bzw. 
25 in die Zylinder (17) als lineare Funktion des Saugrohr- 

A 

drucks (Ps) durch eine Massenbilanz der Luf tmassenstrome 

A A 
• • 

(moK , mzyi) beschrieben wird (Gleichung 2.1) 

- diese ModellgroSen liber eine Dif f erentialgleichung ver- 
knupft werden (Gleichung 2.5), daraus als bestimmende GroSe 

30 zur Ermittlung der tatsachlichen Last der Brennkraft- 

A 

maschine der Saugrohrdruck (Ps) berechnet wird (Gleichung 
2.7) und 

- aus dem linearen Zusammenhang (Gleichung 2.4) zwischen be- 
rechnetem Saugrohrdruck (Ps) und der ModellgroSe fur den 
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Luftmassenstrom (m Zyl )in den bzw. in die Zylinder (17) durch 
Integration die in den bzw. in die Zylinder (17) einstro- 
mende Luftmasse (mzy) erhalten wird. 

5 2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, 

daS das vom Last sensor (12; 13) gemessene Lastsignal 

{m D K_L MM ; Ps_s) zur Korrektur und damit zum Abgleich der Mo- 

A 

dellgrofien (m Zy i) in einem geschlossenen Rgelkreis herangezo- 

10 gen wird, wobei das Lastsignal {mok_uui ; Ps_s) als Fuhrungs- 
groSe des Regelkreises dient. 



15 



3. Verfahren nach Anspruch 2, 

dadurch gekennzeichnet, 
daS der Abgleich im stationaren und/oder instationaren Be- 
trieb der Brennkraf tmaschine durchgefuhrt wird und dabei das 
Ubertragungsverhalten des Lastsensors (12; 13) beriicksichtigt 



wird. 



20 4. Verfahren nach Anspruch 2, 

dadurch gekennzeichnet, 
daS jedetn gemessenen Wert des Drosselklappenof fnungsgrades 
ein Wert eines reuzierten Querschnitts der Drosselklappe 
{Ared) zugeordnet ist und der Abgleich der ModellgroSen durch 

25 Korrektur des reduzierten Querschnitts [Ared) durch eine Kor- 
rekturgroSe {A Ared) derart erfolgt, daE die Regelabweichung 
zwischen FiihrungsgroSe und entsprechender ModellgroSe mini- 
miert wird. 



30 5. Verfahren nach Anspruch 4, 

dadurch gekennzeichnet, 

daS der reduzierte Querschnitt (Ared) aus stationaren Messun- 
gen am Mot orpruf stand ermittelt wird und in einem Kennfeld 
eines Speichers der elektrischen Steuerungseinrichtung abge- 
35 legt ist. 

6. Verfahren nach Anspruch 1, 
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dadurch gekennzeichnet, 
daS bei de£ Darstellung der ModellgroSe filr den Luftmas- 
senstrom (moK) an der Drosselklappe (11) eine in der Durch- 
fluSgleichung (Gleichung 2.2) vorhandene Durchf lufif unktion 
5 (i|/) in einzelne Abschnitte {i « l...k) unterteilt wird und 
diese Abschnitte durch Geradenabschnitte angenahert werden, 
wobei fur die Steigung (mi) und fur das Absolutglied <ni) der 
jeweiligen Geradenabschnitte als Funktion des Verhaltnisses 
von Saugrohrdruck (Ps) und Umgebungsdruck (Pu) bestimmt sind 
10 und in einem Kennfeld abgelegt sind. 

7. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daS die Steigung (y,) und das Absolutglied (y 0 ) der linearen 
15 Funktion fur die ModellgroSe fur den Luf tmassenstrom in den 
bzw. ind die Zylinder (mz y t) abhangig von mindestens einem 
der Parameter, Drehzahl der Brennkraf tmaschine, Zylinderzahl , 
Saugrohrgeometrie, Tempera tur der Luft (T s ) im Ansaugrohr 
(10) und Ventilsteuerzeichen festgelegt sind. 

20 

8. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daS die Parameter duch stationare Messungen am Motorpruf stand 
ermittelt werden und in Kennf elder abgelegt sind. 

25 

9. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daS die in den Zylinder einstromende Luftmasse (mzyt) durch 
die Beziehung 

30 

A A ^ 

mzyi[N-\] + mz y i[N] 

J 



berechnet wird, mit 

Abtastzeit oder Segmentzeit 
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24 

A 

mzyi[N] ModellgroSe des Luf tmassenstromes wahrend des 
aktuellen Abtastschrittes oder Segments 

m Z yi[N-\] ModellgroEe des Luf tmassenstromes wahrend des 
vergangenen Abtastschrittes oder Segmentes. 

5 

10. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 
daB die in den bzw. in die Zylinder einstromende Luftmasse 
imzyi) fur einen, bezuglich der aktuellen Lasterf assung zum 
10 Abtastzeitpunkt [N] in der Zukunft liegenden, bestimmten Pra- 
diktionshorizont (H) geschatzt wird durch Schatzung des ent- 
sprechenden Druckwertes nach folgender Beziehung: 

mz yt [N + H]=T A ^ yi -jij [N]+{H-0.S)-^Ps [N-l]+P s [*]| + y 0 

15 

mit 



T A : Abtastzeit oder Segmentzeit 

H: Pradiktionshorizont, Anzahl der in der Zukunft 

liegenden Abtastschritte 
20 Yi : Steigung der linearen Gleichi^ng 

y 0 : Absolutglied zur Bestimmung von mz y i 

N: aktueller Abtastschritt 



11. Verfahren nach Anspruch 10, 
25 dadurch gekennzeichnet, 

daS die Anzahl (H) von Segment en, fur die das Lastsignal fur 
die Zukunft geschatzt werden soil, drehzahlabhSngig festge- 
legt ist. 
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